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Abstract 

We study in this paper the scenario where the dark matter is constituted by Majo- 
rana particles which couple to a light Standard Model fermion and an extra scalar 
via a Yukawa coupling. In this scenario, the annihilation rate into the light fermions 
with the mediation of the scalar particle is strongly suppressed by the mass of the 
fermion. Nevertheless, the helicity suppression is lifted by the associated emission 
of a gauge boson, yielding annihilation rates which could be large enough to allow 
the indirect detection of the dark matter particles. We perform a general analysis of 
this scenario, calculating the annihilation cross section of the processes XX ~~ ffV 
when the dark matter particle is a SU(2)l singlet or doublet, / is a lepton or 
a quark, and V is a photon, a weak gauge boson or a gluon. We point out that 
the annihilation rate is particularly enhanced when the dark matter particle is 
degenerate in mass to the intermediate scalar particle, which is a scenario barely 
constrained by collider searches of exotic charged or colored particles. Lastly, we 
derive upper limits on the relevant cross sections from the non-observation of an 
excess in the cosmic antiproton-to-proton ratio measured by PAMELA. 
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1 Introduction 



Among the various proposals to characterize the dark matter in our Universe, the scenario 
where the dark matter is constituted by weakly interacting Majorana particles stands 
as the most promising one. In this scenario, thermal scatterings of Standard Model 
particles in the early Universe can produce a relic density of dark matter particles which 
is of the correct order of magnitude, when the interaction strength is of the order of 
the weak interaction strength and the dark matter mass is about 1 TeV. Furthermore, 
this scenario has the appealing feature that the dark matter particle might be directly 
detected in underground detectors, indirectly detected in cosmic ray detectors, gamma- 
ray and neutrino telescopes, and directly produced at the LHC. 

In this paper we will focus on the possibility of indirectly detecting Majorana dark 
matter particles via their self-annihilation in the Milky Way dark matter halo. Concretely 
we will study the annihilation process into two fermions and one gauge boson which, 
under some circumstances, can have a non-negligible or even a larger cross section than 
the lowest order annihilation process into two fermions. More specifically, the s-wave 
contribution to the thermally averaged cross section for the 2 — ?• 2 process is helicity 
suppressed by the mass of the final fermion, while the p-wave contribution is suppressed 
by the small velocity of the dark matter particles in the Milky Way halo. In contrast, for 
the 2 — > 3 process the s-wave contribution is no longer suppressed, due to the associated 
emission of a vector in the final state. As a result, the 2—7-3 processes can even have 
a larger cross section than the 2—7-2 processes as the lifting of the helicity suppression 
can compensate the suppression due to the additional coupling a em /n, provided the 
mediating scalar particles are not too heavy. 

Dark matter annihilations into two fermions and one photon with the mediation of a 
heavy scalar particle were first studied in the framework of the Minimal Supersymmetric 
Standard Model in neutralino annihilations [U[2], and explored in a number of subsequent 
papers [3]. It was also pointed out that this process not only could have a sizable cross- 
section but also produces a gamma-ray with a very peculiar spectral shape which, if 
detected in gamma-ray telescopes, could be unequivocally identified as being originated 
in dark matter annihilations jH [5] . If the mediating scalar particle is electrically charged 
it must necessarily carry hypercharge and therefore the annihilation process must also 
produce weak gauge bosons and, in turn, antiprotons. This process has been studied in [6l 
[7] employing a toy model where the dark matter particle is a singlet under the Standard 
Model gauge group and the mediating particle is a SU{2) L doublet. The constraints on 
the annihilation cross section in this model from the non-observation of an excess in the 
cosmic antiproton-to-proton fraction measured by PAMELA were derived in [7]. Lastly, 
annihilations into quarks, so that the emission of a gluon is allowed in the final state 
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were discussed in [21 [8] . 

In this paper we aim to extend this analysis, considering various toy models where the 
dark matter particle is a SU(2)l singlet or doublet, which couples to the left-handed or 
right-handed leptons or quarks of the first generation. For each case, we will calculate the 
cross sections for the different 2—7-3 processes and we will calculate the constraints on the 
cross sections from the PAMELA measurements of the antiproton-to-proton fraction [9]. 

The paper is organized as follows. In section [2} we discuss in detail the cross sec- 
tions for the various electromagnetic and electroweak internal bremsstrahlung processes 
occurring for SU (2) l singlet as well as doublet dark matter particle coupling to leptons. 
The corresponding cases, when assuming a coupling to quarks, are discussed in section 
[3j In section H] we present the constraints on the cross sections from the PAMELA mea- 
surements of the antiproton-to-proton fraction, and translate them into upper limits on 
an astrophysical boost factor for the various toy models. Finally, we conclude in section 
[5j Our full analytical results for the cross sections of all 2 — > 3 processes considered in 
this work can be found in the Appendix. 



2 Dark matter coupling to leptons 

We consider an extension of the Standard Model by one Majorana fermion, x, which 
we assume to constitute the dominant component of dark matter in the Universe, and 
one scalar particle, 77, which mediates the annihilation process into light fermions. The 
Lagrangian is 

n p 1 f* 1 f* 1 /'fermion , /'scalar /n i\ 

L — *-SM + L, x + U v + Mnt +Mnt • \ ZA > 

Here, £sm is the Standard Model Lagrangian which includes a potential for the Higgs 
doublet $, V = m\&$ + |Ai($ t( l>) 2 . On the other hand C x and £ n are the parts of the 
Lagrangian involving just the Majorana fermion x and the scalar particle 77, respectively, 
and which are given by 

2 2 x (2.2) 

C v = (D^)t(D^) - m\r)S - -A 2 (^) 2 , 

where _D M denotes the covariant derivative. Lastly, £ m ™ lon and CJ^ W denote the fermionic 
and scalar interaction terms of the new particles to the leptons and to the Higgs doublet. 
These terms depend on the details of the model and will be discussed case by case below. 

The quantum numbers of the relevant Standard Model particles under the gauge 
group SU(3) C x SU(2) L x U(1) Y are: e R = (1, 1, -1), L e = (1, 2, -±), $ = (1, 2, \). On 
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the other hand, the quantum numbers of the dark matter particle and the scalar r\ are 
constrained in our setup by the requirement that the dark matter particle is colorless 
and electrically neutral, and by the requirement that a Yukawa coupling to the leptons 
(either left-handed or right-handed) is invariant under the Standard Model gauge group. 
We will assume in the following that the dark matter particle only couples to the first 
generation of leptons, which can be ensured by postulating that the extra scalar particle 
rj carries electron lepton number L e = —1, while the dark matter particle does not carry 
lepton number. Lastly, in order to guarantee the stability of the dark matter particle, we 
impose a Z 2 discrete symmetry under which \ an d V are odd while the Standard Model 
particles are even. 

In the scenarios of interest for this paper, the intermediate scalar is electrically 
charged and could possibly lead to experimental signatures in collider experiments. Pre- 
cise measurements of the invisible decay width of the Z boson at LEP set the upper 
bound Ar inv < 2.0 MeV [10], which rules out the existence of exotic charged scalar 
particles with mass below 40 GeV [11]. Furthermore, the OPAL collaboration searched 
for an excess with respect to the Standard Model expectations of dilepton events with 
missing energy induced by the production of exotic scalar charged particles which de- 
cay into an electron and an invisible particle (in the framework of supersymmetry, the 
production of selectrons which decay into an electron and the lightest neutralino). The 
non-observation of an excess in a sample of 680 pb _1 of e + e~ collisions at center-of-mass 
energy between 192 GeV and 209 GeV, leads to the lower bound m n > 97.5 GeV, assum- 
ing m v — moM > 11 GeV [12]. A similar search was undertaken by the L3 collaboration 
using a sample of 450 pb _1 collisions at ^fs = 183 — 209 GeV, resulting in the lower 
bound > 94.4 GeV assuming m v — tudm > 10 GeV [13], by the ALEPH collaboration 
using a sample of 207 pb _1 collisions at y/s = 204 — 209 GeV, resulting in m, > 95 GeV 
assuming — ttidm > 15 GeV [H] and by the DELPHI collaboration using a sample 
of 609 pb" 1 collisions at ^ = 192 - 208.8 GeV, resulting in m„ > 94 GeV assuming; 
m v — m DM > 15 GeV, and m, > 98 GeV assuming m v — m DM > 5 GeV and m DM < 60 
GeV [15] . For smaller mass splittings the detection efficiency is significantly reduced and 
the lower bounds derived by the LEP experiments can be avoided. 

In the remainder of this section we present a classification of models, according to 
the charge of the dark matter particle under SU(2)l- 

2.1 SU(2)l singlet dark matter 

When the dark matter particle is a SU(2)i singlet, its hypercharge must be zero in order 
to render an electrically neutral particle, hence the gauge quantum numbers must be 
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X = (1, 1, 0). On the other hand, the gauge quantum numbers of the scalar rj depend on 
whether the dark matter particle couples to the right-handed electron singlet or to the 
left-handed electron doublet. 

If the dark matter has a Yukawa coupling to the right-handed electron singlet and a 
scalar field rj, then gauge invariance requires rj = (1, 1, 1). With these quantum numbers 
the only interaction terms in the Lagrangian are: 

4r ion = -fxe R V + hx. , 

In the Minimal Supersymmetric Standard Model (MSSM) this possibility is realized if x 
is the bino and rj is the right-handed selectron, e_R. 

On the other hand, if the dark matter only couples to the left-handed electron doublet 
then rj = (1, 2, — |). The interaction terms are then: 

£ fer m ion = _/x(L e «^) + h.C. = -fx("eLV° ~ ^ + ) + h.C. , 

^t lar = -A a (* t *)(»? t '?) - A 4 ($tr?)(r?t$) . 

In the MSSM, this possibility is realized if x is the bino and rj the left-handed selectron 
doublet, L e . 

After the electroweak symmetry breaking, the mass of the electrically charged scalar 
rj^~ is given by m 2 ± = m\ + \%v 2 EW - If the dark matter couples to the left-handed electron 
doublet, there exists also a neutral scalar with mass given by m 2 = m 2 , + (A 3 + Xi)v EW . 
We will assume that m x < m^i, such that the Majorana fermion x is stable and can 
constitute the dark matter. The interactions lead to a thermal production of x m the 
Early Universe that is compatible with the WMAP value if the coupling is of order one 
/ ~ 0(1) and the masses lie between the weak and the TeV scale [Ml [7]. 

The annihilation of dark matter in the Milky Way today can proceed via the 2 — > 
2 annihilation channels XX ~^ e< 3, as well as XX ~^ y v- However, the cross-sections 
ov2~^2 = a + bv 2 are highly suppressed because the s-wave contribution a oc m 2 is 
helicity suppressed, while the p-wave contribution bv 2 is suppressed by the dark matter 
velocity v ~ 10~ 3 c in the Milky Way halo. The helicity suppression can be lifted by 
emitting an additional spin-1 particle in the final state [H [21 El El El- The lifting of 
the helicity suppression can compensate the suppression due to the additional coupling 
OLem/it provided the mediating scalar particles are not too heavy, typically m v i < 5m x [7J. 
Therefore, the dominant annihilation processes are 2—7-3 channels, like XX ~^ 7 e ^, 
XX -> Zee or XX -> Weu. 

Note that one could similarly consider a coupling of the dark matter particle to the 
leptons of the second or third generation. The cross-sections for the 2 — > 3 processes and 



the constraints from the antiproton flux that will be discussed later are independent of 
the lepton flavor to a good accuracy. Due to the larger masses of the /x and r leptons, 
the helicity suppression of the 2—7-2 annihilation cross-sections is less pronounced than 
for a coupling to electrons, while the velocity suppressed contributions to the 2—7-2 
annihilation cross-sections are flavor-independent. 

The annihilation mode XX ~~ >* leads to a gamma ray signal with a pronounced 
peak at the dark matter mass, that is potentially observable by the Fermi-LAT (see [17] 
for a recent discussion) and by current and future IACTs (4j [5]. On the other hand, 
the annihilation channels involving weak gauge bosons yield a primary contribution to 
the cosmic flux of antiprotons. In the following, we will analyze the relative strength of 
these channels, and in Section 4 we will derive upper limits on the cross-section from the 
PAMELA measurement of the antiproton to proton ratio [9] . 

Coupling to right-handed electrons 

Let us first consider the possibility that the dark matter particle x couples to right- 
handed electrons. Since the mediating particle rj carries hypercharge and electric charge, 
the 2—7-3 annihilation channels 7ee and Zee are possible [6j [7J- The branching ratios 
compared to the 2 — > 2 annihilation cross section are shown in the upper panel of Fig. [TJ 
For dark matter masses 3> far above the Z-threshold, the ratio of the cross 

sections for the electromagnetic and the electroweak bremsstrahlung processes approach 
constant values given by the ratio of the respective coupling constants: 

vv(xx Ze e) ■ ^(xx -> 7ee) = tan 2 (9 w ) ~ 0.30. (2.5) 

The general formulas for the double differential cross sections and for arbitrary dark 
matter masses are given in the Appendix. In order to obtain a gauge-invariant result 
it is important to take into account the diagrams for which the gauge boson is emitted 
off the internal line and off the final state particles (later on, when considering doublet 
dark matter, also contributions from initial state radiation have to be included). The 
corresponding Feynman diagrams are also shown in the Appendix. In the following, 
we will refer to all these processes as internal bremsstrahlung (IB). As noted before, 
these comprise also the contributions from final state radiation in general. Note that 
the contributions from soft and collinear emission, which are in principle logarithmically 
enhanced, are typically negligible in this context because they are helicity suppressed, 
like the 2—7-2 processes. Instead, the dominant contribution arises from the diagrams 
where the gauge boson is emitted either off the internal line, or from a final state particle 
with an off-shell intermediate state (see Ref. [7] for a detailed discussion). 

The dependence of the cross sections on the mass m v is shown in the upper right panel 
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of Fig. [TJ It is apparent from the figure that the branching ratio of the 2 — > 3 processes is 
largest when is close to m DM - Furthermore, for /i = rn^/m'^ )M 3> 1, the cross section 
of the 2 — y 3 processes fall of as l//i 4 , while the 2—7-2 cross sections scale like l//i 2 . The 
former remain dominant as long as the dark matter mass and the mass of the mediating 
particle are of comparable size, roughly < 5m DM - The qualitative properties discussed 
here are common also to most other cases considered below. However, there are some 
quantitative and also qualitative differences which we will stress in the following. 



Coupling to left-handed electrons 

If the dark matter particle x couples to the left-handed electron doublet, it can give rise 
to annihilations into final states involving 7, Z or W bosons. Note that this case has 
been discussed in detail in Ref. [TJ. We will briefly review it here for completeness. The 
branching ratios compared to the 2 — > 2 annihilation cross section are shown in the lower 
part of Fig. [TJ For m DM 3> ^p, and assuming that m v o = m v ±, the asymptotic values 
are again given by the ratios of the appropriate couplings: 



vv(xx -> Zee) 
vv(xX Zvv) 
w(xx -> Wev) 



w(xx^iee) = cot 2 (29 w ) = 0.41, 
av{xx^iee) = = 1.41, (2.6) 

av(xx^jee) = = 4.32. 

Here crv(xx ~~ >* Wev) = crv(xx ~~ >* W~ev) + av(xx ~~ >" W + ev). 

Generically, one expects a non-zero mass splitting of the neutral and charged com- 
ponents of 77 induced by the breaking of the electroweak symmetry, m^„ — m 2 ± = \av\ w . 
Compared to the degenerate limit, the branching ratios get modified due to two effects. 
First, the masses in the t-channel propagators of the mediating particles corresponding 
to the charged and the neutral component of 77 differ from each other. Second, the mass 
splitting opens up a new channel, namely the annihilation into longitudinally polarized 
VT-bosons. In the limit m v i 3> m^M 3> ^f 2 -, the branching ratios are approximately given 
by 

avjxx -> Zvv) _ 1 a4 

vv(xx -> 7ee) sin 2 (29 w ) /4 ' 

W{XX Wev) _ 1 4 [ 1 1 5 m 2 _ J 

crv(xx -> 7*0 " sin 2 (M f [ L + 8 M 2 , l/i± M \ ' ^ <} 

where /i, = rn 2 i /m^ )M and fi = (fi± + fio)/2. The ratio av(xx ~~ >" Zee)/av (xx ~~ >* 7*), 
in contrast, is not affected by the mass splitting. The emission of longitudinal W bosons 
also leads to a spectrum that is harder compared to the case m^o = m v ± [TJ. Analytical 
expressions for the double differential cross sections, from which the spectra can be easily 
obtained, are given in the Appendix. 



SU(2)i singlet DM coupling to e R m 
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Figure 1: Ratio of three-body and two-body annihilation cross-sections, for electro- 
magnetic IB, ov(xX ~^ l e ^)/ av (xX e ^)i anc l f° r the electroweak IB channels 
XX ~~ >" Zee, XX ~~ Zvv and XX ~^ Wev. The latter denotes the sum of W~ev and W + ev. 
The top and bottom rows show the case of SU{2) L singlet dark matter coupling to the 
right-handed electron, and to the left-handed electron doublet, respectively. The left col- 
umn shows the dependence on the dark matter mass for fixed ratio m^±/mDM = 1-5, 
while the right column shows the dependence on the mass of the mediating scalar particle 
r\ for m DM = 300GeV. For the relative dark matter velocity we use v = 10 _3 c. 
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2.2 SU(2) L doublet dark matter 



In this case the dark matter doublet must have one electrically neutral component, which 
is achieved by postulating that the hypercharge is ±1/2. This matter content, however, 
leads to gauge anomalies which can be canceled by introducing another SU{2) L doublet 
with opposite hypercharge. Then, the minimal model with SU(2)l doublet fermionic 
dark matter must contain the new fermions Xi = (1, 2, — |), Xi = (1, 2, |), both charged 
under the Z 2 discrete symmetry. In the case of the MSSM these two particles can be 
identified with the two higgsinos. 

Under these assumptions, the only gauge invariant and Z 2 invariant fermionic mass 
term in the Lagrangian is My\io2Xi + h.c, which generates identical tree level masses 
for Xii Xii X25 xt ■ Quantum corrections induced by the Standard Model gauge bosons 
generate a mass splitting between the charged and the neutral component of the mul- 
tiplet {m x ± — m x o) ra <i ~ 0.34GeV [18], inducing the decay of the former into the lat- 
ter. Therefore, this toy model predicts the existence of two stable particles, candidates 
of dark matter, x± = (Xi i X-D/v^, which will annihilate, among other channels, 
X+X-7 X-X- ~^ e + e~V, with V a vector. Since we are interested in the general fea- 
tures from annihilations of SU(2) doublet dark matter particles and not in constructing 
a fully realistic model, we will assume in what follows that only one of these, x+ 01 X-i 
is present in our Universe today. This can be achieved by postulating a mass splitting 
between them, so that one of them decays into the other at very early times, for instance 
by introducing the dimension five operators 



r /'fcrmion J 

^-^mass ^ 



c 1 (xiia 2 $*)(<f>Ua 2 x c 1 ) + c 2 (x 2 $)($ T X2) + c 3 (X2*)($ f ^2X^ 



h.c. , 

(2.8) 

with A a mass scale larger than the Higgs vacuum expectation value and q coefficients of 
order one. Let us denote the dark matter mass eigenstate as x-, the heavier neutral state 
as x' an d the charged component by x^ 1 - The mass splittings induced by the dimension 
five operator are given by, 

v 2 

Sm± = m x ± - m x = -^(c 3 + |ci - c 2 |) , 

5rriQ = m x i — m x = E ^ \c\ — c 2 | , (2.9) 

up to corrections of order 0(v% w / '(A 2 m x )). Up to the same order, the mass eigen- 
states are related to the two doublet fields by xi = Pl({x' + e x)/V^, X~) an d X2 = 
Pl{{x~Yi (x' ~~ e x)/ A /2), where e = sgn(ci — c 2 ), x an d x' are Majorana fields and x~ 
is a Dirac field. We will assume in the following that the radiative corrections to the 
mass splittings can be neglected compared to the ones induced by the dimension five 
operators. 



By decomposing the gauge interactions of xi an d Xi into mass eigenstates one obtains 

+ exi^A^x' + ecot(26 w )x-rZ^ X - • (2.10) 

Note that this coincides with the interactions of the neutralino in the higgsino limit 
within the MSSM. These interactions give rise to the dark matter annihilation channels 
into a pair of weak bosons with cross-sections given by 

av xx ^ ww = ^ {ml + ml± _ M , wf ^-M w /ml 



3 4 ml - Ml 



av 



a /nr. — m 7 i ~ 

„„ = * Z / 1 ^2/^2 

XX ^ ZZ 64vrc^(m2 -Mf) 2 



Ml/ml . (2.1i; 



As is well-known, these cross-sections can be altered substantially for dark matter masses 
in the TeV range, and if the mass splittings are of order GeV or below, by the multiple 
exchange of weak bosons among the fermions in the initial state, analogous to Sommerfeld 
enhancement in electrodynamics [19]. We will assume here that the mass splittings are 
large enough, such that the effect of Sommerfeld enhancement is in a perturbative regime, 
and comment on its impact below. As we will see this requires A < lOTeV. 

In the present work, we are motivated by the observation that internal 
bremsstrahlung can lift the helicity suppression of fermionic final states. In fact, as we 
will discuss below, the annihilation XX ~~ >* e + e~V can be under certain conditions as im- 
portant as the gauge processes XX ~~ >" WW, ZZ. Let us start by discussing the fermionic 
interactions analogous to the SU (2) singlet case. In particular, the dark matter can cou- 
ple to the right-handed electron singlet and a scalar field 77 = (1, 2, —5)) yielding the two 
following interaction terms in the Lagrangian: 

4nt mi ° n = /(Xi^)eH + h.c. , 
Ct lar = -A 3 ($ t $)(r/ t r / ) - A 4 ($tr / )(r / t$) . 

Alternatively, the dark matter particle can couple to the left-handed electron doublet 
and a scalar field 77 = (1, 1, 1), 

4~ = /(Xi^ 2 ^ + h.c, 

In a supersymmetric context, the scalars can be identified with L e and en, respectively. 
Note that one could in principle consider additional scalar particles that lead to analogous 
couplings involving Xi- We do not consider this possibility in the following. 

The relic abundance produced by the thermal freeze-out is determined by the 2 — ?- 2 
cross- sect ions arising from gauge and Yukawa interactions. If the latter are subdominant, 
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an abundance in accordance with the WMAP value can be achieved for dark matter 
masses of order TeV [18]. By adding the Yukawa interactions, and adjusting the coupling 
/, it is in principle possible to obtain the WMAP value also for even higher dark matter 
masses. However, in the following we will not restrict the range of the dark matter mass 
or the coupling / in order to determine the constraints arising from the measurements 
of the antiproton flux in a way that is independent of the production mechanism. 

We will first discuss the relative size of the electromagnetic and electroweak 
bremsstrahlung in analogy to the case of singlet dark matter, and then the relative 
importance of fermionic and diboson final states. Throughout, we will assume that 

— m x ^> 5m±, Sm , and use the notation moM = Tn x in analogy to the singlet 
case. 



Coupling to right-handed electrons 

If the mediating scalar rj has quantum numbers rj = (1,2, — |), it leads to annihilations 
of dark matter into right-handed electrons. For jtidm 3> - and m v 3> ?tidm, and 
taking only annihilations mediated by the scalar r\ into account, the branching ratio of 
electromagnetic to electroweak bremsstrahlung is given approximately by 

av{ XX ^Zee) : av( XX -»- lee) ~ ^Z^^+^lzl (2 .14) 

where \i = mjL / m^u, and sw and cw are the sine and cosine of the weak mixing angle. 
More general analytical expressions are given in the Appendix. The neutral component 
77 plays no role because it could only lead to the production of right-handed neutrinos, 
which are absent in the SM. Note that the branching ratio increases 

with the mass of the mediating particle to the fourth power, /i 2 oc m^ ± . The reason 
for this behaviour is that the dark matter particle, being an SU(2) 1 doublet, couples also 
to the Z boson. Therefore, annihilation to Zee can occur also via initial state radiation. 
The latter leads to a non-zero contribution to the s-wave cross-section already at the 
l/yU 2 -level, while electromagnetic bremsstrahlung occurs at order l//i 4 . A similar result 
has been obtained within an effective operator approach for Wino-like dark matter in 
[20] . Note that, for very large values of /i, eventually the p-wave contribution to 'jee will 
dominate over the s-wave contribution to 7ee. However, the former also scales like l//x 2 . 
This means the ratio of electroweak and electromagnetic cross sections saturates for very 
large values of /x, which can be estimated roughly as ji ~ (5v 2 ln 2 (mDM/'^e))~ 1 ^ 2 ~ O(30) 
for v = 10 _3 c and wdm ~ 10 2 GeV. However, since a strong gamma signal with spectrum 
peaked at high energies is produced only for [i ~ 0(1), we will not discuss this case in 
further detail. The full dependence of the branching ratios on the dark matter mass and 
the mass of the mediating scalar particle is shown in the upper part of Fig. In the right 

11 



part, it is also shown that the p-wave contribution to the annihilation into 7ee becomes 
important when /i is large. 

Since the electroweak bremsstrahlung is strongly enhanced compared to electro- 
magnetic bremsstrahhmg even for moderate values \i > 1.5, it is important to investi- 
gate whether higher-order contributions can lift the l//i 4 suppression of electromagnetic 
bremsstrahlung. It turns out that this is indeed the case, when considering the correc- 
tions arising from Sommerfeld enhancement. Here, we will estimate the leading effect 
when the enhancement is perturbatively small, following Ref. [211 EH]. In general the 
s-wave amplitude A xx ^sm for annihilation into some SM final state can be written as 

<A XX -+SM = SqA° xx ^ sm + So^x'x'-*-^ S ±^x+x"^SM (2.15) 

where the amplitudes -4° denote the tree-level amplitudes for annihilations of the vari- 
ous components of the doublet, and = d x ipi(x)\ x ^o are enhancement factors. They are 
related to the wave-functions (fi(x) for radially symmetric (s-wave) two-fermion initial 
states, where x = r ■ mvMv/c is a dimensionless variable related to the spatial separation 
r of the fermions. The wave-functions are solutions of a set of coupled Schrodinger equa- 
tions in the presence of a Yukawa potential oc e~ mvr jr that is generated by the exchange 
of vector bosons V — W, Z, 7 among the fermion pair [21], [22]. For the dark matter masses 
we are interested in we can safely apply the low- velocity limit v/c ~ 10~ 3 <C Mw/m-DM 
and assume that S = rriDuv 2 <C 26m±,25mQ. The latter condition implies that the 
charged and heavier neutral components of the doublet cannot be produced on-shell, 
such that their wave-functions decay exponentially at large separations. At leading order 
in the gauge couplings one then finds the approximate solution 

S ~ 1, S Q ~ —= =, S± ~ -/= = . (2.16) 

\J2s\w M z + V2m DM dm 2V2% M w + x/2m DM Sm± 

The approximation can be expected to hold if |s' |, \ s±\ <C 1. For the range of dark matter 
masses we are interested in, this is safely the case if 5m±, Smo > C(l)GeVx (mDM/TeV) 
(see e.g. [IE]). Here, we will assume that this inequality holds and thus that the tree- 
level cross sections yield a reliable estimate. Nevertheless, for the annihilation involving 
electromagnetic bremsstrahlung, the annihilation proceeding via an intermediate charged 
fermion pair, XX ~~ X + X~ ~ >* l e ^i can be important. Concretely, we find that for ^> 

WDM 

vvjxx -> -»• lee) _ , 2 ° v (x + X~ ~j» 7ee)U=o _ . 2 50/i(/i- 1) + 57 

crv{xx -> 7ee) ~~ ± ov(xX -> 7ee) ~ ± 15 

Here the Sommerfeld enhanced cross section is normalized to the leading order cross- 
section, and the cross section for x + X~ ~~ >* 7 e ^ should be evaluated for an initial state 
with total spin zero [21] (note that, for the same reason, the channel XX ~~ X + X~ ~~ ^ 
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Figure 2: Ratio of three-body and two-body annihilation cross-sections, for the case of 
SU(2)l doublet dark matter coupling to the right-handed electron, and to the left-handed 
electron doublet, respectively. The left column shows the dependence on the dark matter 
mass for fixed ratio m„±/fflDM = 1-5, while the right column shows the dependence on the 
mass of the mediating scalar particle rj for moM = 300GeV. Also shown is the correction 
from Sommerfeld enhancement (SE) for mass splittings 5m± = l,10GeV between the 
charged and the lightest neutral components of the dark matter doublets. For the relative 
dark matter velocity we use v = 10~ 3 c. 



ee is helicity suppressed although x is a Dhac particle). Thus, as expected, we find 
that for the annihilation via a charged intermediate state the cross section scales like 
1/yU 2 instead of l//i 4 as for the tree-level contribution. Therefore, if \i is large enough, 
it may compensate the suppression factor |s±| 2 , and yield a significant contribution 
to the annihilation via electromagnetic bremsstrahlung. The influence of Sommerfeld 
enhancement is shown also in Fig. [2] for mass splittings 8m± = 1, lOGeV. Note that 
analogous corrections exist also for electroweak bremsstrahlung. However, since their 
leading order cross sections scale already like l//i 2 , these will be small corrections even 
when /i is large. 
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Figure 3: Cross-sections of various annihilation channels for the case of SU(2)l doublet 
dark matter coupling to the left-handed electron doublet, for m^/moM = 1.5 (top) and 
1.01 (bottom). In both cases, the Yukawa coupling is chosen as / = 2. The annihilation 
cross-section into Wev and Zee has two contributions, one from t-channel exchange of 
the scalar r\ and one from annihilation into WW or ZZ with a subsequent decay of 
one of the gauge bosons. The contribution from the latter are shown also as grey lines. 
Above the threshold for WW / ZZ, the latter processes are taken into account already by 
the 2 — ?- 2 cross-sections, and are shown here for illustration only. For the relative dark 
matter velocity we use v = 10 _3 c. 
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Coupling to left-handed electrons 



If the mediating scalar 77 has quantum numbers 77 = (1, 1, 1), it leads to annihilations of 
dark matter into left-handed electrons. For m^M ^> ^ and 3> m^M, we find the 
following ratios 

av{xx^Zee) : av( XX -+ 7*0 ~ 

°" s h' c w ( 2 18) 

av(xX^Weu) : ot,( xx -> lee) ~ 5 ° M 6 ^ )+63 



Annihilation into Z^P does not occur at tree level, because the quantum numbers of 77 
imply that it couples the dark matter particle only to the charged lepton. Also, since 
77 is a SU(2)l singlet, only transversally polarized W bosons are produced via initial 
as well as final state radiation provided that 5m±,5mo <C M\y- The reason for the /i 2 - 
dependence of the ratios of cross sections is due to initial state radiation, as discussed 
above. The full dependence of the branching ratios on the dark matter mass and the 
mass of the mediating scalar particle is shown in the lower part of Fig. [2j 

When including Sommerfeld corrections, the annihilation channels Zvv and ^vv 
appear. In particular, the latter yields a contribution that scales like l//_i 2 instead of 
l//i 4 as for 7ee. Its cross section is given by 

<rv(xx ->• X + X~ -» 7^) _ \ c± \ 2 ° v ( x+x ~ — ^)l 5 =° ~, \ c± \ 2 ^^ 2 + 12 (2 19) 
<r»{xx -> 7ee) ~~ ± crv(xx -> 7*) ~ ± 15 

The Sommerfeld corrections shown in the lower panel of Fig. [2] refer to the sum of the 
electromagnetic annihilation cross sections uv(xx ~~ >* l^e) + uv{xx ~^ X + X~ ~~ ivv). 

The relative size of the cross sections XX ~~ >" WW, ZZ and XX ~~ >* Vee, Vvv, Wev 
depend on the ratio of the Yukawa coupling / to the gauge coupling g as well as the 
ratio of the mass of the mediating particle 77 and the dark matter mass. In general, due 
to the lifting of helicity suppression, the 2 — > 3 processes are much less suppressed than 
the corresponding 2 — ^ 2 annihilations into fermions. However, the 2 — > 2 annihilations 
mediated by gauge interactions are also not helicity suppressed and therefore generically 
dominate over the 2—7-3 channels. Nevertheless, there can be some cases when the 
latter are important. For example, this can be the case if the dark matter mass is in the 
range M w < ttidm ^ 2M W , such that the diboson states are kinematically disfavored. 
Another possibility is a rather large value for the coupling /. In Fig. [3] the various cross 
sections are shown for / = 2. In this case, the cross sections of the 2 — ^ 3 channels arising 
from electroweak bremsstrahlung are the dominant annihilation channels provided that 
m r)/ m DM i$ 1.01. For even larger /, the latter restriction can be relaxed. We note that 
couplings of that size are required for dark matter masses in the multi-TeV range, when 
imposing the relic density constraint from thermal freeze out. In order to determine 
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constraints from the antiproton flux, we will consider both the case that the annihilation 
channels XX ~~ >* WW, ZZ or XX ~~ >" Vee, Wev dominate. 

3 Dark matter coupling to quarks 

The analysis of the annihilation of dark matter particles into quarks is completely paral- 
lel to the annihilation into leptons discussed in the previous section, the main difference 
being the inclusion of a color quantum number in the final fermion states and in the 
intermediate scalar state. As a result, in addition to the electromagnetic and electroweak 
bremsstrahlung processes, a new annihilation channel arises where a gluon can be ra- 
diated off the final quark states or off the internal colored scalar state. Being both the 
gluon and the photon massless gauge bosons, the resulting spectra will be identical [2]. 
However, the cross section for the gluon internal bremsstrahlung will be enhanced com- 
pared to the electromagnetic internal bremsstrahlung by the larger coupling constant 
and by the color factor, therefore we expect a larger impact of the 2 — ^ 3 processes in 
scenarios where the dark matter couples to quarks compared to scenarios where the dark 
matter couples to leptons, and in particular a larger impact of the present measurements 
of the antiproton-to-proton fraction on the constraints on the couplings of the model. 

Scenarios where the dark matter particle couples to quarks and to a colored scalar 
particle are strongly constrained by experiments searching for exotic colored particles. 
The LEP constraint on the invisible width of the Z boson, Ar inv < 2.0 MeV [10], allows 
to set the absolute lower limits m, > 44 GeV if r\ = (3, 2, |) or (3, 2, |) and m Tj > 33 
GeV if rj = (3,2, — |), corresponding to the searches for left-handed quark doublets, 
right-handed up quarks and right-handed down quarks, respectively |llj . 

Recently, the ATLAS collaboration has reported in [24] the results for the search of 
squark and gluinos using final states with jets and missing transverse momentum using 
1.04 fb _1 of data taken in proton-proton collisions with a/s = 7 TeV at the Large Hadron 
Collider. In this analysis it is considered a simplified supersymmetric model with R-parity 
conserved containing only squarks of the first two generations, q, a gluino octet, g, and 
the lightest neutralino, x?, while all other supersymmetric particles are assumed to be 
very heavy. In this simplified scenario, the supersymmetric particles are produced in pairs 
and the final states produce more than 2, 3 or 4 jets plus missing energy. Concretely, the 
production of two squarks qq is followed by the decay q —> qx\, thus producing at least 
two jets plus missing energy. In contrast the production of two gluinos gg is followed by 
the decay g — > qqXi, which yields in the final state at least four jets plus missing energy. 
Lastly, the associated production of one squark and one gluino qg yields at least three jets 
plus missing energy. The non observation of an excess in any of these channels above the 
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Standard Model background can be translated into constraints on the (m g ,mq) plane, 
assuming m v o = 0. 

X-l 

Our toy model for dark matter corresponds to the simplified SUSY model considered 
by the ATLAS collaboration in the limit where the gluino is also very heavy and is not 
kinematically accessible to the LHC with yfs = 7 TeV. Hence, we conclude that the 
non-observation of an excess over the Standard Model background of dijet events with 
missing energy in the present LHC data translates into a lower bound on the colored 
scalar state of our toy model > 875 GeV at 95% c.l. It is important to note that this 
lower bound assumes that the mass splitting between \ an d t] is large enough to produce 
jets passing the requirements for the transverse momenta (p? > 130 GeV for the leading 
jet and pt > 40 GeV for the second jet). Then, this stringent lower bound on the colored 
scalar mass can be avoided if x an d r\ present a degenerate mass spectrum, as discussed 
in [25], concretely when m v — ?tidm < 130 GeV, so that the dijet event does not pass all 
the cuts required by the ATLAS analysis. 

Searches for colored scalar states were also undertaken at the Tevatron and at LEP. 
The searches at the Tevatron by the CDF and DO [27] collaborations employ similar 
cuts as the ATLAS analysis described above, and have by now been superseded. On the 
other hand, the searches at LEP, despite limited by the smaller center of mass energy 
and by the smaller luminosity, employed a smaller cut for the jet transverse momentum 
and are relevant for our analysis. Concretely, the L3 collaboration has presented limits 
on the squark masses searching for an excess in dijet events with missing energy in e + e~ 
collisions at center of mass energies between 192 GeV and 209 GeV with an integrated 
luminosity of 450.5 pb _1 [13J. The non-observation of an excess with respect to the 
expected Standard Model background translates in our toy model into the lower bound 
m v > 97 GeV for m v — m DM > 10 GeV. 

As a summary, we conclude that the toy model with a Majorana dark matter particle 
which couples to the quarks and a colored scalar, t], via a Yukawa coupling is in agreement 
with the present searches of new physics if 

• m v > 875 GeV for any 7Bdm- 

• 97 GeV < m v < 875 GeV, if m v - m DM < 130 GeV. 

• 33 -44 GeV < m v < 97 GeV, if m v - m DM < 10 GeV. 

Following the same scheme as in the case of annihilations into leptons, we will analyze 
the features of various dark matter scenarios coupling to quarks according to the charge 
of the dark matter particle under SU(2)l- 
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3.1 SU(2)l singlet dark matter 



This choice of the SU(2)l charge requires that the gauge quantum numbers of the dark 
matter particle must be x — (1> 1)0) in order to render an electrically neutral particle, 
while the quantum numbers of the intermediate colored scalar rj depend on whether the 
dark matter particle couples to the right-handed up quarks, Ur = (3,1,1), the right- 
handed down quarks, dn = (3, 1, — ~), or to the quark doublet, qi = \ d L ) = (3, 2, |). 

The fermionic interaction term when the dark matter couples to the right-handed 
up quark reads: 

4;™ = -fXURV + h.c. , (3.1) 

which requires quantum numbers for the intermediate scalar particle rj = (3, 1, — |), while 
the scalar interaction term is given by Eq. fl2.3p . Similarly, a Yukawa interaction between 
the dark matter and the down quark requires rj = (3, 1, |). In the MSSM these particles 
correspond to the right-handed up and down squarks, respectively. 

On the other hand, the fermionic interaction Lagrangian of the dark matter particle 
to the quark doublet reads: 

4n r t mi ° n = -fx(q L i*2r)) + h.c = -fxM - d L rj d ) + h.c. , (3.2) 

with rj = (Zd) = (3,2, — |), which in the MSSM corresponds to the left-handed squark 
doublet. The scalar interaction Lagrangian is given in Eq. fl2.4p . 

In this case, the helicity suppression of the 2—7-2 annihilation channels XX ~^ 
urUr, dndn, QlQl: can be lifted by the associated emission of photons, weak gauge 
bosons or gluons together with the light quarks. Let us analyze the relative strength of 
these channels for each of the scenarios. 

Coupling to right-handed up quarks 

In this case the particle mediating the dark matter annihilations carries hypercharge, 
electric and color charge, therefore the annihilation channels '-jurUr, ZurUr, qurUr are 
possible. We show in Fig. 0], upper plot, the corresponding cross sections relative to the 
helicity suppressed cross section for XX ~^ u rMr- 111 the left plot we present the ratio of 
cross sections for dark matter masses between 50 GeV and 5 TeV and different values 
of the mass splitting — m DM = 100 GeV, 50 GeV and 10 GeV, to study the impact 
of the constraints from collider searches of exotic colored particles. It is apparent from 
the plot that the smaller the mass splitting, the larger is the relative cross section of the 
2—7-3 processes, especially for light dark matter particles. On the other hand, in the 
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limit ttidm ^ ~2~ the ratios are fairly independent of the mass splitting and take the 
values: 

<™{XX QUrUr) : <rv(xx -»• JUrUr) = 3a s (m DM ) / a ern ~ 38.4, 

<rv(XX -* ZurUr) : <rv(xx -»• JUrUr) = t&n 2 (6 w ) = 0.30. (3.3) 

Exact formulae for the various cross sections can be found in the Appendix. For the 
numerical values given here and below, we have evaluated the strong coupling constant 
m D M = 300GeV for illustration. 

Moreover, we show in the upper-right plot the dependence of the ratio of cross 
sections on the mass of the intermediate colored scalar particle for a fixed dark matter 
mass todm = 300 GeV. The maximal value of the ratio of cross sections, which is as large 
as 4.6 x 10 4 for XX ~~ >* 9 u rUr, is reached when ~ tjdm- We note that this is precisely 
the region of the parameter space which is most difficult to constrain at colliders, as the 
jet produced in the decay of the colored scalar particle is too soft to be triggered. As in 
the leptonic case, for fi = m^/mp M 3> 1, the cross sections of the 2 — )• 3 processes scale 
as 1/yU 4 , while that of the 2 — y 2 process scale as l//x 2 . 

Coupling to right-handed down quarks 

The results for the annihilations into right-handed down quarks are completely analogous 
to the results for the annihilations into right-handed up quarks presented above, the only 
difference being the different hypercharge (and electric charge) of the intermediate scalar. 
As a consequence, the cross sections for the annihilations XX ~^ {Z,l)dRdR are a factor 
of 1/4 smaller than the corresponding cross sections for XX ~^ i)urUr. In particular, 
when mDM 3> ^ the cross sections for annihilations into right-handed down quarks 
satisfy the relations: 

<rv(xx ->■ gdRd R ) : (tv(xX -> ld R d R ) = 12a s (m DM )/a em = 154, 

w(xx -> Zd R d R ) : crv(xx -> ld R d R ) = t&n 2 (9 w ) = 0.30 . (3.4) 

Coupling to left-handed quarks 

When the dark matter particle couples to the left-handed quarks, a new annihilation 
channel is open, involving W bosons which can be radiated off the internal colored scalar 
or off the final fermions legs. The ratios av (xx ~ ^ QlQl) / '&v(xx ~^ QlQl) are shown in 
Fig. m lower-left plot, in the limiting case mjj = m v ± for dark matter masses between 
50 GeV and 5 TeV and different mass splittings between the dark matter mass and the 
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intermediate scalar mass, m v — ttidm = 10, 50 and 100 GeV. Compared to the coupling 
into right-handed quarks, it is noticeable the enhancement of the branching ratio into 
electroweak gauge bosons, arising from the additional channel involving W bosons and 
both up and down quarks, as well as a stronger coupling of left-handed compared to 
right-handed quarks to the Z-boson. On the other hand, the branching ratio into gluons 
is identical to the right-handed case, while the branching ratio for ^qhQh differs by a 
factor of 8/5 (2/5) with respect to ^urUr (jdndji) due to the different electric charges 
of the particles involved and due to the doubling of diagrams. In the limit m^u 3> ^f- 
the cross sections satisfy the relations: 



ov{xx -> gqq) 

w(xX -> Zqq) 
vv{.XX -> Wqq') 



vv(xx ->■ IQO) = 24a s (m DM )/(5a em ) = 61.4, 

<rv{xx^iw) = ^W^r 1 = 3 - 02 ' (3 - 5) 

°v{xx^iqq) = 9/(54,) = 7 - 79 , 

with av(xx ~^ Wqq') = av(xx ~^ W + (Ilu l ) + crv(xx ~~ W^u^d^). 

The dependence of the branching ratios with the mass of the intermediate scalar 77 
is shown in Fig. HJ lower-right plot. 

In a more realistic scenario, the two weak isospin components of the intermediate 
scalar particle will not be degenerate in mass, but will have a mass splitting proportional 
to the order parameter of the electroweak symmetry breaking: m^ — m^ ± = \±v\^. 
As already discussed for dark matter particles coupling to leptons, in this situation the 
t-channel propagators in the annihilation diagrams have poles at different masses and, 
more importantly, the new channel with annihilations into longitudinally polarized W- 
bosons opens up, resulting in an enhancement of the branching ratio for XX ~~ WqiqL- 
In the limit m v i ^> m^M ^ we ^- n( ^ ^ na ^ ^he cross sections for the 2—7-3 annihilation 
cross sections satisfy: 

°v(xx ^ gw) _ 12 a s {m DM ) /i« 4 + /^ 4 



VV{XX ->■ iqq) "em 4/i~ 4 + /i d 4 ' 

av(xX -> Zqq) = (3 - As 2 w f^ + (3 - 2s 2 w ) 2 ^ A 
<™{XX ->• iqq) sin 2 (2^)(4/x- 4 + ^ 4 ) 

av(xX^Wqq') 9((/i u + fJL d )/2)-'- 



(3.6) 



5 m E)M / n5 



™(xx -»• jqq) s 2 w (4fi- 4 + n/) 

where \X{ = m^jm^y^. Analytical expressions for the differential cross sections can be 
found in the Appendix. 



3.2 SU(2) L doublet dark matter 

The case of doublet dark matter coupling to quarks can be discussed in close analogy 
to the case with a coupling to leptons. In particular, the gauge quantum numbers of the 
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Figure 4: Ratio of three-body and two-body annihilation cross-sections. The top and 
bottom rows show the case of SU{2)l singlet dark matter coupling to the right-handed 
up-quark, and to the left-handed quark doublet, respectively. The left column shows the 
dependence on the dark matter mass for fixed mass splitting m v — ttidm = 10, 50 and 
lOOGeV, while the right column shows the dependence on the mass of the mediating 
scalar particle rj for m^M = 300GeV. For the relative dark matter velocity we use v = 
10~ 3 c. 
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Figure 5: Ratio of three-body and two-body annihilation cross-sections. The top and 
bottom rows show the case of SU(2)i doublet dark matter coupling to right-handed up- 
and down-quarks with equals strength (/ = /'), and to the left-handed up quark, respec- 
tively. The left column shows the dependence on the dark matter mass for fixed mass 
splitting m v — m^M = 10, 50 and lOOGeV, while the right column shows the dependence 
on the mass of the mediating scalar particle r\ for 7«dm = 300GeV. For the relative dark 
matter velocity we use v = 10 _3 c. 
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Figure 6: Cross-sections of various annihilation channels for the case of SU(2)l doublet 
dark matter coupling to the left-handed up quark, for m,±/moM = 1.5 (top) and 1.01 
(bottom). In both cases, the Yukawa coupling is chosen as / = 2. For the relative dark 
matter velocity we use v = 10~ 3 c. 
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dark matter doublets Xi,2, their mass eigenstates and gauge interactions are identical. 
Differences arise only with respect to the coupling to colored scalar particles 77. We 
will again discuss the different possibilities. The dark matter particle can couple to the 
right-handed quarks via a mediating scalar 77 = (3, 2, 1/6), 

4n r t mi ° n = f{xii°W*)d R + f'(x2ta 2 V*)u R + h.c. . (3.7) 

In contrast to the leptons, two coupling terms are allowed by the symmetries, coupling 
the dark matter to the up- and down quarks, respectively. Alternatively, the dark matter 
particle can couple to the left-handed quark doublet. In this case, two choices for the 
quantum numbers of the mediating scalar field are possible. For rj = (3,1,-1/3), the 
interaction term reads 

4T ion = fiXi^Dv + h.c. = f(-j=& + tx)dl ~ X^l) + h.c. , (3.8) 

which couples the dark matter particle x to the left-handed down quarks. The second 
possibility is to choose 77 = (3, 1, 2/3). Then the interaction term reads 

4~ = f(x^ 2 q c L ) V + h.c. = /( - -Ltf - ex)ul + x±d|) + h.c. , (3.9) 

and we obtain a coupling of dark matter to left-handed up quarks. In a supersymmetric 
context, the scalars can be identified with qi, cIr, and ur, respectively. These interactions 
will lead to annihilation into light quarks via internal bremsstrahlung of either a photon, a 
weak gauge boson, or a gluon. For the photon and the gluon, the corresponding diagrams 
are identical to those obtained for singlet dark matter, and therefore the cross sections 
are also the same. However, electroweak bremsstrahlung can proceed also via initial 
state radiation. As in the leptonic case, this leads to a considerable enhancement of the 
annihilation into Wqq' and Zqq. We will now discuss the branching ratios for each case. 



Coupling to right-handed quarks 



The annihilation of dark matter into right-handed quarks can be mediated by a coloured 
scalar rj with quantum numbers (3,2,1/6). There are two possible couplings, / and 
/', which correspond to a coupling to ur and to cIr, respectively. Here, we assume for 
simplicity that both components of rj are degenerate in mass. The branching ratios in 



the limit m v » m DM » ^ are given by 



w(xx ->■ 79?) 
vv(xx -> 79?) 
vv(xx -> 79?) 



vv(xx -> Wqq') 



(4/ 4 +/' 4 ) 



12(/ 4 +/' 4 ) « s (-m DM ) 
4/ 4 +/' 4 



f 4 F(4s 2 w /3)+f' 4 F(-2s 2 w /S) 



<i0swc 2 w 
Hff') 2 50 M 2 +12 
4/ 4 +/' 4 15^ ' 



(3.10) 
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where F(x) = 50/i(/i + x) + 15(1 — x) 2 — 3, fi = m^/m 2 ^^ and the cross sections denote the 
sum of annihilation processes into up and down type quarks. As in the leptonic case, the 
ratio of cross sections of electroweak to electromagnetic processes depends quadratically 
on fi because of the different scaling of the cross sections with \i due to initial state 
radiation. In the upper right part of Fig. [5] the dependence of the cross sections on the 
dark matter mass is shown for three different mass splittings m,- ttidm = 10, 50, lOOGeV, 
assuming that / = /'. Note that the annihilation into weak bosons can be as strong or 
even stronger than the annihilation into gluons. The reason is that the former can proceed 
also via initial state radiation, which leads to a parametric enhancement compared to the 
annihilation into gluons for large values of fi. It turns out that even for moderate values of 
\x the additional channels due to initial state radiation enhance the electroweak processes 
significantly. Note, however, that the branching ratio into W bosons gets suppressed when 
the ratio of couplings \ f / f'\ deviates from unity. The dependence of the branching ratios 
on is shown in the upper left part of Fig. |5j As expected, the cross sections for gqq and 
jqq fall off as l//i 4 oc 1/m®, until the p-wave contribution dominates for m v > 4mDM- 
The electroweak processes scale with 1/fi 2 , because of the contribution from initial state 
radiation. Since also the 2 — 2 cross section scales like l//i 2 , the ratio approaches a 
constant. In Fig. [5j also the leading effect of Sommerfeld enhancement is shown, for the 
same choice of parameters as discussed in the leptonic case. As for a coupling to leptons, 
the main effect is to enhance the branching fraction into 'jqq due to the annihilation 
process XX ~^ X + X~ ~^ 1QQ- The corresponding cross section can be inferred from the 
formulae given in the Appendix. 



Coupling to left-handed quarks 

There are two possibilities, rj can have the quantum numbers (3, 1,2/3) or (3, 1, —1/3). 
In the first case it mediates a coupling of dark matter to Ul, and in the second to d^. 
Lets consider both possibilities separately. For rj = (3,1,2/3), the branching ratios in 
the limit m v 3> m DM ^> read 

<?v{xx^guu) : av(x X juu) ~ 3a s (m DM ) / a em , 

av(xx^Zuu) : <jv{ XX ^ juu) ~ 3F(1 - As 2 w / 3) / (80 s 2 w c 2 w ) , (3.11) 

av( X X^Wqq') : av( X X ~> iuu) ~ 3(50/i(/i - 1) + 63)/(80s^) . 

For t] = (3, 1, —1/3), on the other hand, the branching ratios are given by 

w(xX^gdd) : av(xx^ldd) — 12a s (m DM )/a em , 

av(xx^Zdd) : <jv( X X^jdd) ^ 3F(2s 2 w /3 - l)/(20s 2 w c 2 w ) , (3.12) 
av(xx^Wqq') : av( X X -> idd) ~ 3(50/i(/i - 1) + 63)/(20 S 2 w ) , 
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where the function F(x) is the same as defined above. The dependence on the dark 
matter mass and on the mass of the mediating scalar is shown in the lower part of Fig. [5] 
for the case t] = (3, 1,2/3). The behaviour is qualitatively similar to the one discussed 
above. In addition, the absolute values of the various cross sections are shown in Fig. [6] 
for the same choice of parameters than for the case of a coupling to leptons. Due to the 
color enhancement and the annihilation channel into gluons, the annihilation processes 
mediated by the scalars rj are important over the whole dark matter mass range for both 
choices m,/moM = 1.5 and 1.01. 



4 Antiproton-to-proton ratio and observational con- 
straints 

Two dark matter particles at the position f can annihilate producing antiprotons at a 
rate per unit of kinetic energy and volume given by: 

CM^EH^, (4.D 

Z '"DM f al 

where (av) f is the thermally averaged cross-section multiplied by the velocity in the 
annihilation channel /, p(r) is the distribution of dark matter particles in the Milky 
Way, where r denotes the position of the dark matter particle with respect to the center 
of our Galaxy, and dNp/dT is the energy spectrum of antiprotons produced in that 
channel per unit of kinetic energy. We will assume for simplicity a spherically symmetric 
distribution, and calculate the antiproton flux assuming a radial dependence given by 
either the Isothermal, NFW or Einasto profile with parameters specified in [7]. The 
spectrum of antiprotons is obtained using the event generator PYTHIA 8.1 [28] interfaced 
with CALCHEP [29l [30]. 

After being produced at the position r, antiprotons propagate through the Milky 
Way in a complicated way before reaching the Earth. Following [21] , we will describe an- 
tiproton propagation by means of a stationary two- zone diffusion model with cylindrical 
boundary conditions. Under this approximation, the number density of antiprotons per 
unit kinetic energy, fp(T,r,t), approximately satisfies the following transport equation: 

° = it = V ' ( ^ (T ' f) V/p_) " V ' &(?)fp) - 2hS(z)T aam f p + Q(T, r) . (4.2) 

The boundary conditions require the solution fp(T,r,t) to vanish at the boundary of 
the diffusion zone, which is approximated by a cylinder with half- height L — 1 — 15 kpc 
and radius R = 20 kpc. The diffusion in the Galactic magnetic field and the convection 
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Figure 7: Upper bounds on the electroweak IB cross-sections av(2 — > 3) obtained from 
the PAMELA data [9] of the cosmic antiproton-to-proton ratio. The plots correspond to 
the constraints on the individual electroweak IB processes XX ~~ >* Weu (where Weu = 
W~ev + W + eu, upper), and XX ~~ ^ 9 U ^ (lower) at 95%C.L. Dashed lines correspond to 
the MIN, solid to MED and dotted to MAX propagation models. The bounds obtained 
assuming an isothermal dark matter profile are shown in red, NFW in black, and Einasto 
in blue. 
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Model 


5 


K (kpc 2 /Myr) 


L (kpc) 


V c (km/s) 


MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.70 


0.0112 


4 


12 


MAX 


0.46 


0.0765 


15 


5 



Table 1: Astrophysical parameters compatible with the B/C ratio that yield the minimal 
(MIN), median (MED) and maximal (MAX) flux of antiprotons. 

term, which accounts for the drift of charged particles away from the disk induced by 
the Milky Way's Galactic wind, are described by the parameterization K{T) = Kq (3 1Z 6 
and V c (r) = V c sign(z) k. The third term accounts for antimatter annihilation with rate 
r ann , when it interacts with ordinary matter in the Galactic disk. In order to take the 
uncertainties related to propagation into account, we will use three sets of parameters, 
compatible with the cosmic boron to carbon flux ratio [32], corresponding to minimum, 
medium and maximum antiproton flux as given in Table [H The flux at the position 
of the solar system is given by $^ S (T) = -£-fp(T, r Q ). Finally, we take the effect of 
solar modulation into account using the force field approximation [331 G2! with solar 
modulation parameter <j)p = 500 MV for our numerical analysis. In order to obtain 
constraints on the dark matter annihilation cross-section, we compute the antiproton- 
to-proton ratio p/p = ($^ lg + $p kg )/$ p using the proton flux of [35]. The background 
flux arises from antiprotons that are produced by spallation of cosmic ray nuclei, mainly 
protons and Helium, on the interstellar medium. We use the backgound flux calculated in 
Ref. [36j based on the two- zone diffusion model, taking into account p-p, He-p, p-He and 
He- He nuclear reactions. The main uncertainties arise from the diffusion parameters and 
the nuclear cross-sections, and are estimated to be in the range of 10 — 25% depending 
on the energy. In contrast, the uncertainty stemming from the knowledge of the primary 
flux of cosmic nuclei and the composition of the interstellar medium are found to be 
subdominant. Note that the prediction for the secondary antiproton flux is consistent 
with the PAMELA data [9]. In order to obtain a conservative exclusion bound we adopt 
the minimal value for the antiproton background as discussed in [36]. Upper limits on the 
cross-sections of the individual annihilation channels, as well as on astrophysical boost 
factors, are then obtained from the PAMELA p/p data [9] using a x 2 -test at 95%C.L. 

We show in Fig. [7] our results for the upper limits on the cross-sections of individual 
IB processes. In order to take the astrophysical uncertainties into account, we compute 
limits for the three propagation models and three halo profiles discussed above. The 
limits on the individual channels do not depend significantly on the mass spectrum as 
long as m^/mj)M ~ C*(i)- For the MED propagation model the limits lie in the range 
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BF max from p/p for singlet DM coupling to L e 




m DM [GeV] 



BF max from p/p for singlet DM coupling to 




rriDM [GeV] 

Figure 8: Upper bounds on the astrophysical boost factor for singlet dark matter coupling 
to leptons. Here we use the MED propagation parameters and the NFW profile. The 
contours correspond to the maximally allowed boost factor log 1Q (BF) from the PAMELA 
p/p data [9]. The coupling / is fixed by requiring that the thermal relic density matches 
the WMAP value. The upper plot shows the case of annihilations into left-handed leptons. 
The solid lines correspond to m^ — m^ ± = v 2 EW and the dashed lines to m v o = m v ± . In 
the lower plot the constraints for dark matter annihilating to right-handed electrons are 
shown. 
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BF maa; from p/p for singlet DM coupling to ur 




m DM [GeV] 

Figure 9: Upper bounds on the astrophysical boost factor for singlet dark matter coupling 
to quarks, obtained under the same assumptions as described in Fig. El The upper plot 
shows the case of annihilations into left-handed up and down quarks. The solid lines 
correspond to — m^ d = v 2 EW and the dashed lines to m Vu = m Vd . In the lower plot 
the constraints for dark matter annihilating to right-handed up quarks are shown. The 
relevant IB process is XX 9QQ- 
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Figure 10: Upper bounds on the astrophysical boost factor for doublet dark matter 
coupling to gauge bosons and leptons. The upper plot shows the case of dark matter 
annihilating exclusively into gauge bosons. In the lower plot we compare the bounds 
that arise from including leptonic final states with a coupling / = 1 and / = 2 to the 
case without leptons. For a better comparison the astrophysical parameters have been 
fixed to a NFW profile and MED propagation in the lower plot. 
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av {xX ~~ ^ Weu) < 10~ 25 cm 3 / sec for ttidm = lOOGeVand 10~ 24 cm 3 / secformDM = ITeV 
for electroweak IB processes involving leptons. For dark matter coupling to quarks the 
IB of gluons is the dominant source of antiprotons. The limit lies in the range av (xx ~~ ^ 
guu) < 3 x 10~ 26 cm 3 / sec for m D M = lOOGeV and 5 x 10~ 25 cm 3 / sec for m D M = ITeV. 
These limits can be compared to the ones recently derived in [17] on the annihilation cross 
section of singlet dark matter particles into right-handed fermions and a photon. The 
non-observation by the Fermi-LAT of an excess of gamma-rays translates into the limits 
av {,XX ~^ 7/ i+ y u ~) ^ 6 x 10~ 28 cm 3 / sec, av (xx ~~ ^ 7T + r~) <6x 10~ 28 cm 3 / sec, av (xx ~~ ^ 
^ybb) <4x 10~ 28 cm 3 / sec for m^M = lOOGeV. For dark matter coupling to leptons, the 
upper limit on the cross section on the process XX ~~ >* Wei/ from the PAMELA antiproton 
data roughly translates, using Eq. 12.61 into av(xx ~^ 7^ + ^~) ^ 2x 10 _26 cm 3 / sec, which 
is weaker than the limit obtained in [TTj from the Fermi-LAT data. For larger masses, 
where there are no limits from gamma-rays, the PAMELA data provide the strongest 
constraints on the model. On the other hand, for dark matter coupling to quarks, one 
obtains from the PAMELA antiproton data and Eqs. (I3.3[ I3.4[ 13.51) the upper limits 
av (xX ~~ >" IQQ) ^ 2 x 10~ 28 cm 3 / sec, 7 x 10~ 28 cm 3 / sec, 5 x 10 _28 cm 3 / sec for couplings 
to down-type quark singlets, up-type quark singlets and quark doublets, respectively, and 
wdm = lOOGeV. Hence, in this range of dark matter masses the antiproton constraints on 
the model can be, depending on the propagation model, competitive with the gamma-ray 
constraints. 

The limits on the cross-sections can be translated into limits on the coupling strength 
of dark matter to Standard Model particles within the scenarios discussed in the previous 
sections. In turn, these can be used to compute limits on an astrophysical boost factor. 
For the case of singlet dark matter, these limits are computed by fixing the coupling 
/ such that thermal freeze out produces a relic density in agreement with the WMAP 
value Q x h 2 ~ 0.11. The required value of the coupling parameter / depends on the mass 
spectrum and can be found in Ref. [T6l 17]. For the case of doublet dark matter, the gauge 
coupling is fixed and therefore the above prescription would be too restrictive. Therefore, 
we do not impose the relic density constraint in that case. Instead, we compute the 
antiproton flux arising from the annihilations of the dark matter particles in the galaxy 
using the physical value of the gauge coupling constant and fixed values / = 0, 1,2 of 
the Yukawa coupling. A limit on the boost factor can then be obtained by requiring that 
the antiproton flux obtained from the annihilations of doublet dark matter for a given 
set of masses and coupling constants, and multiplied by the boost factor, is consistent 
with the PAMELA p/p data at the 95%C.L. Note that the maximally allowed boost 
factor determined in this way can be formally less than unity, which means that the 
corresponding model parameters can be excluded, (old: Instead, we compute the boost 
factor for the case of doublet dark matter by keeping the ratio of Yukawa- and gauge 
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coupling at a fixed value, and rescale all the fluxes according to a common boost factor.) 

The results for singlet dark matter coupling to leptons are shown in Fig. [HJ in de- 
pendence on the dark matter mass and the mass splitting between the scalar particle 
rj mediating the annihilation and the dark matter particle. In Fig. EE the corresponding 
constraints for singlet dark matter coupling to quarks are shown. For leptons, the rele- 
vant channels are XX ~~ >* Weu, XX ~~ >* Zee an d XX ~^ Zvv. For quarks, the antiproton 
flux is dominantly produced by the process XX ~~ >* As expected, the constraints are 
significantly stronger than for leptons, mainly because the internal bremsstrahlung of 
gluons has a larger cross section and the colored particles in the final state lead to a 
more efficient antiproton production. 

In Fig. [10] the results for doublet dark matter coupling to leptons are shown. In this 
case the channels XX ~~ >* W + W~ and XX ~^ ZZ are relevant for most couplings and mass 
splittings. Only for a coupling / > 1 and a mass splitting m,/mDM ~ 1.01 non-negligible 
contributions arise from the leptonic channels XX ~~ >" Weu and XX Zee. It should be 
noted that doublet dark matter with masses ttidm ^ 200GeV can be excluded due to 
the antiproton constraints even without any additional antiprotons produced in leptonic 
channels. This reach can be extended through the inclusion of leptonic channels. As was 
to be expected the constraints worsen significantly for dark matter masses m x > ITeV, 
so that doublet dark matter is virtually unconstrained by antiprotons in the TeV range. 

In general, the antiproton constraints are most stringent for small mass splittings 
m v — ttidu, because the internal bremsstrahlung cross sections are strongly enhanced. 
Note that in this parameter region the bounds from collider searches are weakest, as 
discussed before. 

5 Conclusions 

We have studied the annihilation process of two Majorana dark matter particles XX ~^ 
ffV, with / a light Standard Model fermion and V a gauge boson, which is relevant for 
their indirect detection. Under very general conditions, this process dominates over the 
lowest order annihilation process XX ~~ //> due to the lifting of the helicity supression 
in the s-wave contribution to the cross section of the 2 — )• 2 process by the associated 
emission of a spin 1 particle. 

To keep the analysis as general as possible, we have performed an extensive analysis 
of possible dark matter scenarios, focusing on the plausible case where the dark matter 
particle is a SU(2)l singlet or doublet. We have classifed all scenarios where the dark 
matter particle couples to the first generation of Standard Model fermions and an extra 
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scalar particle, which mediates the annihilations XX ~~ >* f f an d XX ~~ >" ffV. The gauge 
invariance of the Yukawa coupling requires the intermediate scalar particle to carry 
hypercharge and electric charge. Hence, the helicity suppression of the 2—7-2 process 
can be lifted through the associated emission of a photon or a weak gauge boson off 
the final charged fermions or off the intermediate charged scalar, which is more efficient 
when the mass of the intermediate scalar is close to the mass of the dark matter particle. 
In the case of the SU(2)l doublet, the weak gauge boson can also be emitted off the 
initial dark matter particle, resulting in an enhancement of the cross section with respect 
to the SU{2)l singlet case. Moreover, when the two weak isospin components of the 
SU(2)l doublets have a sizable mass splitting, the annihilation cross section into weak 
gauge bosons is further enhanced by the radiation of longitudinally polarized H^-bosons. 
We have provided analytic expressions for the cross sections for all these processes, 
complementing results already existing in the literature, and we have studied numerically 
the relative strength of each of the annihilation channels. If the dark matter particle 
couples to a light quark via a Yukawa coupling with a colored scalar, then not only the 
radiation of a photon or a weak gauge boson can lift the helicity suppression, but also 
the radiation of a gluon. 

Models with exotic charged or colored particles are strongly constrainted by the 
negative searches at accelerators of an excess over the Standard Model expectations of 
dilepton or dijet events with missing energy. We have remarked that the choices of pa- 
rameters where the 2—7-3 processes are most important, namely when the intermediate 
scalar mass is close to the dark matter mass, evade the stringent lower bounds on the 
masses of the exotic particles from accelerator searches, since the lepton or the jet pro- 
duced in the decay of the exotic charged or colored particle is too soft to pass all the 
cuts required by the current analyses. 

Lastly, we have calculated the constraints on the various dark matter scenarios from 
the non-observation of an excess in the cosmic antiproton-to-proton fraction measured by 
the PAMELA collaboration. We have presented these constraints as upper limits on the 
cross section of the relevant annihilation processes, and translated them into constraints 
on an astrophysical boost factor for each scenario and a broad range of dark matter 
masses as well as mass splittings between dark matter and the scalar particle mediating 
the annihilation. 
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Note added 

During the last stages of this work we became aware of another group discussing the 
annihilation of Majorana dark matter particles into quarks and a gluon [37]. 



A Cross sections 

A.l SU(2) L singlet dark matter 

The differential cross-sections for the two-to-three dark matter annihilation processes 
XX ~> V f f m the limit v, mf — > 0, for the case of S , L r (2)^-singlet Majorana dark matter 
X coupling to the SM fermions / via a mediating scalar rjf are given by 



(A.l) 



vdajxx -> iff) C lff a em f\l - x)[x 2 - 2x(l - y) + 2(1 - yf] 
dE^dEf 8ir 2 m^ M (l - 2y - /i f ) 2 (3 - 2x - 2y + /i f ) 2 

vdajxx Zff) C Zf ja em f A 

dE z dE f 87r 2 m£ M (l ~ 22/ - ^/) 2 (3 -2x-2y + ^) 2 ' 

x |(1 - x)[x 2 - 2x(l -y) + 2(1 - yf] 

+ xl[x 2 + 2y 2 + 2xy - Ay}/ A - 4/8} , (A.2) 

vdajxx -> Wff) C Wf f,a em f 4 

dE w dE f 8vr 2 m4 M (l - 2y - /i / ) 2 (3 - 2x - 2y + ^) 2 

x {(1 - x)[x 2 - 2x{\ -y) + 2(1 - y) 2 + 2(2 - x - 2y)Afj] 

+ x 2 [x 2 + 2y 2 + 2xy - Ay + 2(2 - x- 2y)A/x + A^ 2 ]/4 - x^/8 

+ A/i 2 [(l - 2z)/2 - (1 - y)(l -x- y)/(2x 2 )}} , (A.3) 

vda( XX C 9// a s (m DM )/ 4 (l - x)[x 2 - 2x(l - y) + 2(1 - y) 2 ] 



dE^dEf 8n 2 ml M (l - 2y - fi f ) 2 (3 - 2x - 2y + ///) 5 



(A.4) 



Here, x = E v /m DM for 7 = 7, W, y = E f /m DM , x = M v /m DM , /j, f = m 2 f /m^ M , 

Hfi = m 2 f jTnriM> an d = 2(A i /' ~ A*/)- The pre-factors are given by the following 
expressions 

C lf f C zf J C Wff> Cggg 



XX -»• Vf R f R 
XX -> V7l/l 



g 2 iV c g 2 iV c tan 2 (M iV c C F 

_2 at- (t 3 /~g/ sin 2 (6W)) 2 , . 7V e at r 
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where qj and £3/ are the electric charge and the weak isospin, respectively, with q e = — 1 
andt 3e = —1/2. For quarks, one has iV c = 3, and Cp = (Nl—1)/{2N C ) = 4/3. The spectra 
of the vector bosons can be obtained by integrating the differential cross-section over the 
fermion energy, with integration limits given by j^ n } m l max = mDM _ (Ey±\/Ey — M v )/2. 
The total cross-section can be obtained by integrating over the remaining energy with 
limits Ey m = My and Ey ax = muM + M v /(4muM)- The relevant Feynman diagrams 
are shown in part (a) and (b) of Fig. [TT] for the exemplary case XX ~^ Zee. 

For comparison, we also quote the leading contribution to the two-to-two cross section 
for v — y 0, in the limit mf = 0, 

fV 1 + A 

ov{ X X // = J 2 7T— % , (A.5) 
where v is the relative velocity. 



A.2 SU(2) L doublet dark matter 

Here we report the 2—7-3 cross sections for the case of 5 l f/(2)i-doublet Majorana dark 
matter x that arise from a coupling to the SM fermions / via a mediating scalar rjf. 
There can be additional contributions due to 2 — >■ 2 annihilations into WW or ZZ, 
with a subsequent decay of one of the gauge bosons. We do not include them here for 
simplicity. Their size depends on the ratio g/f of gauge interactions and the interactions 
with the scalar rjf. 

The main difference compared to the case of SU(2) l singlet dark matter are addi- 
tional contributions from initial state radiation. For annihilation into 7// or gff, there 
are no such contributions because the dark matter is electrically neutral and uncoloured. 
Therefore, the cross-section is identical to the one from above. The cross-section for the 
annihilation to Zff is given by a sum of ten diagrams: four with initial and final state ra- 
diation, respectively, and two where the Z is emitted from the mediating particle 77. The 
Feynman diagrams for the case XX ~^ Zee are shown in Fig. [TTJ For the cross-section, 
we find in the limit m/, v — > 

vdajxx -+ Zff) N c a em f 4 

dE z dE f 2ix 2 sin 2 (2^)m 4 M (2x - x 2 ) 2 (l - 2y - /i / ) 2 (3 - 2x - 2y + ^ } ) 2 

x {(1 + x 2 jA - x)[x 2 - 2x(l -y) + 2(1 - yf - x 2 /2}C f (x) 2 

+ x 2 (l -y- x/2) 2 [(1 + x 2 /4 - x)(AC f (x) + x 2 ) 

- 4/2-2(1 -y)(l-x-y)]} , (A.6) 

where C f (x) = 1 + fj, f + (x - xl/2)(g{ ± g v ) - Xq/2, and g v = t 3f - 2q f sm 2 (9 w ) and 
g^ — tsf are the couplings to the Z boson. The plus and minus sign applies to annihilation 
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Figure 11: Feynman diagrams contributing to the annihilation channel XX ~^ Zee. For 
a singlet dark matter particle x the diagrams corresponding to final state radiation (a) 
and virtual internal bremsstrahlung (b) contribute. Note that only their sum is gauge 
invariant. For a doublet dark matter particle Xi a l so the diagrams (c) where the Z-boson 
is emitted from the initial state have to be taken into account in addition to (a) and 
(b). Here xi an d X2 refer to the two neutral mass eigenstates as explained in section [2721 
Their mass splitting is assumed to be negligibly small in Eq. (1A.6|) . 
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into left-handed or right-handed fermions, respectively. For the annihilation into Z/l/l 
the corresponding mediating particle rjf has to have the SM gauge quantum numbers of 
fn, and vice versa. 

Similarly, for the annihilation into W bosons, we find 

vdajxx -)• WJjJD N c a em f 

dE w dE f 64vr 2 s\n 2 {9 w )m% u {2x - x 2 ) 2 (l - 2y - n f ) 2 

x |4(1 - x)[x 2 - 2x(l -y) + 2(1 - y) 2 } + 2(1 - x - y)x% 
+ xl[hx 2 - 2x - 2 + 8(1 - y)(l - x))] + 4/4} , (A.7) 
vda(xx -> Wf R f' R ) 2 vda(xx -> Zff) 

A0\ 



dE w dE f w dE z dE f 



M z ^M w J^(ff')2,C f (x)^l+^ f ~x2/2 

(Ai 



The former process can proceed via a mediating particle with the quantum numbers of 
the right-handed partners of either / or of /' and incorporates contributions from initial 
and final state radiation. The latter process is mediated by the doublet (r)f,r)'f) with 
quantum numbers of the left-handed fermion doublet (/, /'), and involves contributions 
from initial state radiation and from diagrams where the W is emitted off the internal 
line. Here we assumed m r]f = m Vft for simplicity. Note that annihilation into W fnf' R is 
only possible for quarks in the SM. It can also exist for leptons if neutrinos are Dirac 
particles. 

In connection to Sommerfeld corrections also annihilations of the charged components 
of the doublets containing the dark matter particle are relevant. Here it is important to 
project out only the contributions where the initial particles have total spin zero. We 
find that 

vda{x + X~ ~± lff)\s=o _ , 2 2 vdajxx ~± Zff) , 

dE^dE f " ASwCw dEwdEf \Mz^,c f ^i +fi+qfX ■ (A.Jj 

For annihilation into right-handed fermions /r, the mediating particle rjf needs to have 
the same quantum numbers as For annihilation into left-handed fermions fi, a me- 
diating particle with quantum numbers of the SU(2)l partner f' L of fi is required. 
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